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ABSTRACT 

We revisit the relation between magnetic-field strength {B) and gas density {p) for 
contracting interstellar clouds and fragments (or, cores), which is central in observa- 
tionally determining the dynamical importance of magnetic fields in cloud evolution 
and star formation. Recently, it has been claimed that a relation B oc is statis¬ 
tically preferred over B oc in molecular clouds, when magnetic field detections 
and nondetections from Zeeman observations are combined. This finding has unique 
observational implications on cloud and core geometry: The relation B oc p^/^ can 
only be realized under spherical contraction. However, no indication of spherical ge¬ 
ometry can be found for the objects used in the original statistical analysis of the B — p 
relation. We trace the origin of the inconsistency to simplifying assumptions in the 
statistical model used to arrive at the B oc p^/^ conclusion and to an underestimate 
of observational uncertainties in the determination of cloud and core densities. We 
show that, when these restrictive assumptions are relaxed, B oc p^/^ is the preferred 
relation for the (self-gravitating) molecular-cloud data, as theoretically predicted four 
decades ago. 

Key words: diffusion - ISM: magnetic fields - ISM: clouds - MHD - stars: formation 
- methods: statistical 


1 INTRODUCTION 

Whether interstellar magnetic fields play a role in the for¬ 
mation of clouds and stars or are affected by cloud and/or 
star formation are old questions that predate the discovery 
of molecular clouds. Mestel & Spitzer (1956) were concerned 
that the estimated very high electrical conductivity of the 
interstellar gas implies that the magnetic flux is frozen in the 
matter and, in the then prevailing picture of star formation 
(Hoyle’s spherical collapse and hierarchical fragmentation), 
magnetic fields would prevent fragmentation and star for¬ 
mation. For this reason, they suggested that ambipolar dif¬ 
fusion (the motion of electrons and ions together with mag¬ 
netic flux relative to the neutrals) would set in at some stage 
and allow a cloud to reduce its magnetic flux and thereby 
fragment and collapse to form stars. 

The first testable prediction of magnetic-field strengths 
in clouds destined to form stars was given by Mestel (1965). 
A cloud collapsing spherically and isotropically while con¬ 
serving its mass (M) and magnetic flux {^b) implies a re¬ 
lation between the magnetic-field strength and gas density: 
M oc pR^ = const, oc BR^ = const'. Hence, eliminating 
the cloud radius R in favor of the density p from the two 


conservation laws yields a scaling between the field strength 
B and the density: B (x . 

Verschuur (1969) summarized the results of Zeeman 
measurements of the field strength in nine Hi clouds 
on a log(R) - log(n) plot, where n is the number density 
(particles/cm^) and concluded, as did other workers, that the 
measurements were in agreement with that theoretical prediction 
- without, however, performing an actual fit to those data. 

Mouschovias (1976a, b) provided the first self-consistent for¬ 
mulation and solution of the problem of the equilibrium of self- 
gravitating, isothermal, magnetic clouds embedded in a hot and 
tenuous intercloud medium. He also considered the implication 
of the contraction of such clouds (or fragments) on the B - p 
relation. In the deep interior of eaeh cloud, he found that the 
ratio of magnetic and gas pressures, a = B"^ISttP, tends to re¬ 
main constant during contraction. In fact, it tends to retain a 
value near unity. For isothermal contraction, P = pnC^, where 
C = (A;BT/m)^/^ is the isothermal speed of sound at tempera¬ 
ture T and mean mass per particle m; the quantity k-B is Boltz¬ 
mann’s constant. Hence, Mouschovias’ result a 1, specialized to 
isothermal contraction, yields B oc For an equation of state 

P oc pO, the relation between B and p becomes B oc po/"^ (see 
review by Mouschovias 1991b). An analysis by Crutcher (1999) of 
a larger sample of clouds than that used by Verschuur, with mea¬ 
sured magnetic-field strengths and number densities (greater than 
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Figure 1. Different geometries of contracting clouds and magnetic fields examined in azu Black arrows represent the direction of the 
magnetic field and bold red arrows the direction of contraction. The B oc relation is uniquely associated with spherical contraction 
and, therefore, has unique observational implications. 


100 cm“^), found a best-fit exponent of 0.47 for the B — p scaling, 
in agreement with the theoretical prediction. Detailed numerical 
simulations by Fiedler & Mouschovias (1992, 1993) predicted a 
slope hv = 0.47 for contracting cores formed by gravitationally- 
driven ambipolar diffusion and evolving from initially magneti¬ 
cally subcritical to supercritical states. More recently, Li et al. 
(2015) inferred n = 0.41 from observations in the massive star 
forming region NGC 6334. For lower densities, although the the¬ 
oretical prediction for an evolving cloud was k 0 (because self¬ 
gravity is not strong enough to compress the cloud perpendicular 
to the field lines), the observational picture has generally been 
less clear, with various scaling exponents derived empirically (see 
Vallee 1997 and references therein; Marchwinski 2012) for sets of 
clouds observed at the same time. 

Crutcher et al. (2010, hereinafter referred to as CWHFTIO) 
revisited the scaling between B and p for a yet larger sample 
of both low-density (primarily HI) and high-density (primarily 
molecular) clouds. They used a Bayesian statistical analysis that 
allowed them to treat nondetections and varying angles between 
the magnetic field and the line of sight, and they optimized a fam¬ 
ily of models consisting of a uniform distribution of magnetic-field 
values between some minimum and maximum, with the maxi¬ 
mum having two distinct branches in its behavior (on a logB - 
logp plot): a flat part at low densities {B independent of p), and 
a power-law scaling at higher densities {B oc p'^), with the expo¬ 
nent Av, the break density po, and the width of the magnetic-field 
strength distribution being the free parameters of their model. 
Their conclusion was that the data prefer k ^ 2/3 and reject 
K ^ 1/2. They took this result to be an indication of “isotropic 
contraction of gas too weakly magnetized for the magnetic field 
to affect the morphology of the collapse.” 

In this work, we examine more closely the observational im¬ 
plications of different geometries of contraction on the B - p rela¬ 
tion. The distribution of forces in a cloud determines its evolution, 
including its geometric shape and the associated B — p relation. 
Although a given (or observed) B - p relation does not necessarily 
imply a unique geometric shape of a cloud, it is nevertheless the 
case that a given (or observed) B — p relation can only be found 
in a very restricted set of geometric shapes, which in turn re¬ 
strict the kind of motions capable of producing those shapes and 
the B - p relation. Here, we test whether the observed shapes of 
the objects (clouds and cores) on which the latest B - p relation 
study (that of CWHFTIO) has been based are consistent with the 


underlying geometries in which the claimed scaling {B oc p^/^) 
could develop. 

In l|2.1l we summarize the B - p relations implied by different 
cloud geometries that could be established by the evolution of 
molecular clouds with frozen-in magnetic fields (no significant 
ambipolar diffusion). Density maps of clouds and cores used in 
CWHFTIO are examined in ^2.21 testing for consistency between 
geometry and the exponent ac. The value k = 2/3 claimed by 
CWHFTIO cannot be reconciled with the observed cloud shapes. 
The source of the discrepancy lies in various assumptions of the 
CWHFTIO analysis, as we show in JS] Relaxing the problematic 
assumptions, we reconcile the observed shapes and the B — p 
relation in El and we show that the value /c = 1/2 is preferred by 
the data over the value n = 2/3. We summarize the conclusions 
in 3S1 


2 CLOUD GEOMETRY AND THE B - p 
RELATION 

In this section we address the connection between the slope of 
the B - p relation and the cloud geometry. First, in ini we 
investigate theoretically the B — p relation implied by different 
geometries of clouds and magnetic fields. Then, in in we exam¬ 
ine the shapes of objects in the CWHFTIO sample and whether 
they are consistent with the claimed slope k, = 2/3. 

2.1 B — p Relations Implied by Different 
Geometries 

2.1.1 Disklike or Slab Cloud with B in the Plane of the 
Disk 

We first consider an oblate (disklike) cloud of half-thickness Zq 
and arbitrarily large radius Rq, uniform density po, threaded by 
a uniform magnetic field Bq in the plane of the disk (see Fig. 
la). A slab-shaped cloud is a special case of this. Contraction 
perpendicular to the plane of the disk to a new half-thickness Z 
increases the density and the magnetic field by the same factor, 
ZqIZ\ hence, B oc p. One should note that, for this kind of con¬ 
traction, the gravitational force per unit mass perpendicular to 
the plane of the disk on a fluid element initially at 2 ; depends only 
on the column density pz, which does not change upon contrac¬ 
tion. However, the magnetic-pressure force per unit mass on that 
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same fluid element, | — VS^/87r|/p, increases upon contraction as 
. Consequently, even if such contraction sets in for some 
reason, magnetic forces will eventually stop it. 

We now consider the same oblate cloud, threaded by the 
same magnetic held, but now the cloud is allowed to contract 
only along the held lines (see Fig. lb); i.e., the half-thickness (or 
polar radius) Zq does not change, but the local extent of the 
cloud along held lines (oc Rq cosO, where 0 is the angle between 
the held lines and a line from the cloud’s centre to the point of 
interest on the rim of the cloud) decreases such that the density 
increases uniformly in the cloud model. The ultimately resulting 
shape is in general one of a prolate, triaxial object, a “filament” 
perpendicular to the held lines. The motions that created this fil¬ 
amentary cloud do not by themselves change the strength of the 
magnetic held; hence, B oc , i.e., B is independent of p. How¬ 
ever, the increased density in the filamentary structure implies a 
stronger gravitational held along the filament, toward the centre 
of the original oblate cloud. Will this filament fragment along its 
length to form at least one more-or-less spherieal eore? 

If this filament contracts as a whole along its length, an 
argument similar to the one in Mouschovias (1976b) shows that 
the magnetic-tension force near the ends of the filament increases 
more rapidly than the gravitational force, so such contraction 
(perpendicular to the held lines) can more easily take place in the 
central part of the filament. If a fragment (or core) is to separate 
out and contract gravitationally in this region, its mass-to-flux 
ratio must exceed the critical value for collapse. 



(Mouschovias & Spitzer 1976). The quantity a is the column den¬ 
sity along held lines (in g cm“^), and B is the magnetic held 
strength. (The constant on the right-hand side of eq. [1] has a 
slight dependence on the geometry of the cloud.) If the original 
size of the cloud, both parallel and perpendicular to the held lines, 
is very large, then the resulting filament will also be very long, 
and it is possible for several fragments to separate out along its 
length, provided that criterion (1) is satisfied for each. The ther¬ 
mal critical mass per unit length of a filament, 2C^ jG (Ostriker 
1964), is not a relevant quantity for the fragmentation of a fila¬ 
ment threaded by a magnetic held perpendicular to its long di¬ 
mension. For each fragment, B oc for as long as the magnetic 
held remains frozen in the matter. Detailed numerical simulations 
(Fiedler & Mouschovias 1992, 1993) showed that ambipolar diffu¬ 
sion sets in in the interiors of initially subcritical molecular clouds 
and leads to an increase in the mass-to-flux ratio toward its criti¬ 
cal value for collapse (see eq. [1]). Prior to establishment of critical 
conditions, the magnetic-field strength increases by at most 30% 
while the density can increase by a very large factor; hence, k, ^ 0. 
After the mass-to-flux ratio reaches its critical value, contraction 
accelerates and proceeds with balance of forces along held lines 
and as rapidly as magnetic forces allow perpendicular to the held 
lines. These are the sufficient conditions for establishment of the 
relation B oc p'^ with n = 1/2. (Actually, the numerical sim¬ 
ulations show that K, = 0.47, meaning that flux-freezing is not 
perfect; some ambipolar diffusion takes place even during the dy¬ 
namical stage of contraction - see Fiedler &: Mouschovias 1993, 
Fig. 9c.) 

2.1.2 Disklike or Slab Cloud with B Perpendicular to the 
Plane of the Disk 

We now consider the same disklike cloud, threaded by the same 
magnetic held as above, but now the held is perpendicular to the 
plane of the disk (see Fig. Ic), with held lines “fanning out” out¬ 
side the cloud, acquiring a characteristic hour-glass shape (not 
shown in Fig. Ic because it is not essential for the present pur¬ 
poses). Such oblate clouds are unavoidable if they form out of 


lower-density interstellar gas in which the magnetic held has an 
ordered component and, locally, the magnetic force is nonnegli- 
gible relative to the gravitational force. The B oc p^D relation 
implied by the gravitational, isothermal contraction of such a 
cloud (both along and perpendicular to the held lines) as well 
as its physical origin were described above in the Introduction 
and at the end of the preceding subsection and they need not be 
repeated here. Fragments can separate out in the interior of the 
cloud and contract dynamically if their mass-to-flux ratio exceeds 
the critical value for collapse given by equation (1). A magnet¬ 
ically subcritical cloud can reach critical conditions because of 
gravitationally-driven ambipolar diffusion, whose modern under¬ 
standing is that it redistributes mass in the central flux tubes 
of molecular clouds, where the degree of ionization is relatively 
small (< 10“^), but it does not lead to flux loss by a cloud as a 
whole (Mouschovias 1979). 


2.1.3 Cylindrical or Filamentary Cloud with B Along the 
Cylinder 

Cylindrical model clouds threaded by a magnetic held along 
the axis of symmetry of the cylinder (see Fig. Id) were stud¬ 
ied exhaustively by Mouschovias & Morton (1991, 1992a, 1992b). 
During lateral contraction (perpendicular to the symmetry axis 
and the magnetic held), both the instantaneous density p and 
magnetic-field strength H of a fluid element at instantaneous 
cylindrical polar radius r increase as 1/r^; hence, B oc p. The 
gravitational force per unit mass at the position of the fluid el¬ 
ement increases as 1/r, while the magnetic force per unit mass 
increases as 1/r^. Consequently, the magnetic forces will stop 
such contraction at some stage. There is neither a critical mass 
per unit length nor a critical mass-to-flux ratio per unit length for 
such self-gravitating filaments to suffer collapse (see Mouschovias 
& Morton 1991, discussion following eq. [38]). The evolution of 
such model clouds was followed numerically by Mouschovias & 
Morton (1992a, b); the density quickly acquires a spatial profile 
approximated by 1/r^. 


2 . 1.4 Cylindrical or Filamentary Cloud with B 
Perpendicular to the Cylinder 

The case in which the magnetic held is perpendicular to the axis of 
the cylinder (Fig. le) is similar to the one discussed in §2.1.1 above 
- the object that started as an oblate (disklike) cloud threaded 
by a magnetic held in the plane of the disk and then contracted 
primarily along the held lines to acquire a triaxial, prolate shape 
perpendicular to the held lines. 


2.1.5 Cylindrical or Filamentary Cloud with B at an 
Angle with respect to the Cylinder 

If the magnetic held is at an angle with respect to the axis of the 
cylindrical cloud, one might think that its component along the 
axis would ensure that there is no critical mass or mass-to-flux 
ratio for lateral collapse, while its component perpendicular to 
the axis would bring in the Mouschovias & Spitzer (1976) critical 
mass-to-flux ratio for collapse and/or fragmentation of the cloud 
in both directions. Then the longitudinal component of B would 
tend to increase as p, while the lateral component would tend to 
increase as Unfortunately, the effect of the magnetic held 

on the evolution of a cloud cannot be deduced correctly by con¬ 
sidering separately the effect of each of its components and then 
superimposing the two results. To visualize the behavior of a pro¬ 
late cloud threaded by a magnetic held at an angle with respect 



4 Tritsis et al. 


to its length, we consider two cases: (a) B is almost perpendicu¬ 
lar to the axis of the cylinder (see Fig. If), and (b) B is almost 
parallel to the axis (see Fig. Ig). 

In case (a), motions along the field lines can, in principle, 
take place and make the extent of the cloud along B as small as 
thermal-pressure and/or turbulent-pressure forces allow. If only 
this evolution took place, B would be independent of p. How¬ 
ever, for self-gravitating fragments (or cores) to form, destined 
to collapse and form stars, contraction perpendicular to the field 
lines has to take place as well. This brings in the critical mass- 
to-flux ratio given by equation (I), and the relation B oc is 
established again. 

In case (b), motions along field lines can transform the cylin¬ 
der into an oblate (disklike) object or even break it up into sev¬ 
eral such oblate objects, with hardly affecting the magnetic-field 
strength. However, as we have already seen, for an oblate object 
to collapse its mass-to-flux ratio has to exceed the critical value. 
Evolution beyond that state results in the relation B oc 
How the fragmentation of an initially critical magnetic flux tube 
affects the B - p relation depends on the number of fragments 
that may form along the flux tube, and is discussed in detail in 
Mouschovias (I99Ic; see § 2.4 and Fig. I therein). 



2.1.6 Spherical Cloud 

The B - p relation implied by spherical, isotropic contraction (Fig. 
Ih) was discussed in JT] The relation B oc p^/^ is unique among 
B - p relations in that only spherical contraction can cause it. If 
contraction along field lines is more rapid than perpendicular to 
the field lines, the exponent n in the relation B oc p'^ becomes 
less than 2/3. If the opposite is true, the exponent n becomes 
greater than 2/3. Thus the recent claim by CWHFTIO that Zee- 
man and density observations, taken in aggregate, yield B oc p^/^ 
has unique, observationally testable implications on the shapes of 
the observed objects (clouds or cores): they must be spherical. 

The severe constraint on the geometry that can produce the 
B oc p^C relation does not get relaxed by a random compo¬ 
nent of B dominating its ordered ^r mean) component inside 
the observed object (cloud or corejJ • While a completely ran¬ 
dom component of the magnetic field does not introduce a spatial 
anisotropy that would tend to destroy the spherical geometry and 
the B oc p‘^C relation, it does not enter the B - p relation, and 
it also does not contribute to the observed line-of-sight B. The 
B - p relation refers to the mean H, and Zeeman observations 
measure that field’s component along the line of sight. 

2.2 Observational Determination of CWHFTIO 
Cloud Shapes 

The conclusion that the exponent n = 2/3 requires spherical 
clouds and cores is difficult to reconcile with the CWHFTIO re- 

^ In any case, at core scales (smaller than about 0.1 pc), the mag¬ 
netic field is both observed (e.g., Girart et al. 2006, Chapman et 
al. 2013) and theoretically expected (Mouschovias 1987a) to be 
dominated by an ordered component. The reason is that, because 
of magnetically-driven ambipolar diffusion, a random component 
of B cannot be sustained on scales smaller than the Alfven length- 
scale, Aa = 7^VA'^'ni^ where va = Bj[Aicp)^I is the Alfven speed 
in the neutrals, and is the slowing-down time of a neutral 
particle due to collisions with ions; for typical molecular clouds, 
Aa = 0.1 pc (Mouschovias 1987a; 1991a). The timescale for 
straightening out field lines that are tangled on a scale I is propor¬ 
tional to C and, for typical molecular cloud densities (10^ cm“^), 
magnetic fields (30/rC), and core sizes (0.1 pc), is much smaller 
than the free-fall time, (e.g., see Mouschovias et al. 2011). 


Figure 2. Archival Herschel-SPIUE 250pm dust continuum emis¬ 
sion map of the core NGC2024. The star represents the coordi¬ 
nates of the core as reported by Falgarone et al (2008). The inset 
image is the outcome of our data processing and the black lines 
represent the principal axes of the core. 


suit at high densities. Although the shapes of molecular clouds 
and their fragments are an important part of the debate on the 
process that regulates cloud evolution and star formation, spher¬ 
ical clouds and cores are not a contender in either the theoretical 
or the observational arguments. 

According to theory, if ambipolar diffusion is the agent 
mainly responsible for core formation, then the cores are ex¬ 
pected to be flattened along the magnetic field because magnetic 
forces act only perpendicular to the field lines. Thus, the expected 
shapes are oblate (Mouschovias 1976b; Fiedler & Mouschovias 
1993), although not necessarily axisymmetric (Basu & Ciolek 
2004; Ciolek & Basu 2006; Kudoh & Basu 2011). If cores form as 
the result of converging turbulent flows (e.g., see review by Mac 
Low & Klessen 2004), then they are expected to have random, 
triaxial shapes, with a slight preference for prolateness (Gammie 
et al. 2003; Li et al. 2004). Predominantly spherical cores and 
clouds are not expected by any formation mechanism. 

Observationally, the issue is how to distinguish between 
oblate and prolate intrinsic shapes from two-dimensional pro¬ 
jections. Spherical clouds and cores are straightforward to spot 
even in projection: their two-dimensional projections are circles 
and the aspect ratio of the projected shapes is very sharply 
peaked at 1. Early work on core shapes, assuming axial symmetry, 
seemed to favour prolate cores (Myers et al. 1991; Ryden 1996). 
However, subsequent investigations which relaxed the axisymme- 
try assumption have consistently yielded triaxial, preferentially 
oblate core shapes (e.g., Jones, Basu & Dubinski 2001; Jones & 
Basu 2002; Goodwin, Ward-Thompson & Whitworth 2002), in¬ 
dependently of tracer or core sample. Tassis (2007) and Tassis 
et al. (2009) found strong indications for triaxial, preferentially 
oblate cores. Spherical cores and clouds on the other hand are not 
common in nature. 

However, it may still be possible that the molecular sam¬ 
ple used in the GWHFTIO study does contain an unusually high 
fraction of spherical objects, which would be consistent with their 
finding of H oc To test for this possibility, we have calcu- 
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Figure 3. 250yL4m dust continuum emission maps from the Herschel Science Archive. Black lines represent the principal axes of the cores. 
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Figure 4. 500/i,m dust continuum emission maps from the Herschel Science Archive. Black lines represent the principal axes of the cores. 
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IRAM04191 aspect ratio= 0.55 L1521-E aspect ratio= 0.76 LI 524-4 aspect ratio= 0.55 
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Figure 5. line emission maps of the J=l-0 transition from the FCRAO survey. Black lines represent the principal axes of the 

cores. 


lated aspect ratios of plane-of-the-sky maps for as many of the 
objects of CWHFTIO as we could find appropriate data. 

We used dust continuum emission images from the online 
data archive of the Herschel Space Telescope (Pilbratt et al. 2010) 
observed with SPIRE at 250 /rm and 500 fim. to create maps of 
18 of these objects, and line emission data of the J=l- 

0 transition from the FCRAO survey of the Taurus molecular 
cloud (Goldsmith et al. 2008) for another 9 of them. The maps 
were centered on the coordinates given by Troland & Crutcher 
(2008) and Falgarone et al. (2008) as referenced by CWHFTIO. 
The angular sizes of the maps were 3 to 4 times the typical sizes 
of cores. The method we used to calculate the aspect ratios is 
based on the first and second moment of the fiux density and 
is described in detail by Tassis et al. (2009). In certain cases, 
cores, as projected on the plane of the sky, are closely spaced; 
for this reason, fiux from a nearby core can severely affect the 
aspect ratio of the object under examination. To avoid this, we 
performed a visual inspection to select a region that contains 
only one core. Finally, in order to remove the background which 
can also affect the analysis, we set a threshold and calculated the 
aspect ratio only from the pixels whose intensity was greater than 
the mean intensity of the final region (e.g., see Fig. [2]and inset). 
The resulting maps of the cores, along with their principal axes 
and their aspect ratios, are presented in Figures [3]-[5] 

Even visual inspection of these maps reveals that these ob¬ 
jects do not appear to be spherical. Indeed, only 4 of these ob¬ 


jects have an aspect ratio consistent with a spherical geometry. 
The mean value of the aspect ratios as computed here is 0.63, 
suggesting that in our sample of cores the preferred geometry is 
a fiattened, oblate one (Fig. [6]). We therefore conclude that the 
shapes of the cores in the CWHFTIO sample are not consistent 
with the spherical geometry implied by the B oc relation. 

In our theoretical analysis of shapes, we considered only the 
case of a pure dataset following a specific geometry and evolution¬ 
ary path. It is also conceivable that a mixture of different object 
geometries, each with its own B — p relation, could yield a value 
of Tv different from the values characterizing the individual ob¬ 
jects, including possibly tc = 2/3. In order for such a scenario to 
be realized, objects with both smaller values of k, (i.e., 1/2 or 0) 
and greater values of tc (i.e., I) need to be present in the sample. 
As discussed in gni values of n greater than 2/3 can be pro¬ 
duced by one-dimensional contraction of a disklike (oblate) cloud 
with its magnetic field parallel to the plane of the disk, or the 
lateral contraction of a long cylindrical cloud with its magnetic 
field along the cylinder (see § 2). Both of these cases were found 
to be unlikely in a combined study of cloud shapes and magnetic 
fields using data for 32 clouds surveyed by the Hertz polarimeter 
(Tassis et al. 2009, see their Figs. 2b and 2e). The latter possi¬ 
bility is also contradicted by our study of the CWHFTIO core 
shapes alone, as it would require a different aspect ratio distribu¬ 
tion than the one shown in Figure[6] The CWHFTIO aspect ratio 
distribution is similar in shape to the one of the sample studied 
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Figure 6. Distribution of the aspect ratios of the cores for which 
we could find observational data. The distribution peaks away 
from unity, in contradiction to what is expected from spherical 
geometry. 


by Tassis (2007), and which was shown to be strongly peaked 
around intrinsically oblate cores. Even a uniform shape distribu¬ 
tion (equal number of oblate and prolate objects) would require 
a much greater fraction of aspect ratios between 0 and 0.4. It 
is therefore unlikely, given what we know regarding core shapes 
and the orientation of their magnetic fields, that the CWHFTIO 
sample contains enough cores that evolve either as B oc or 
as 5 oc p to effectively “pull” the observed logB - logp relation 
away from a slope of 0.5. 


3 TRACING THE SOURCE OF THE 
DISCREPANCY BETWEEN 
NONSPHERICAL SHAPES AND 
CWHFTIO’S = 2/3 

We have found in JJJthat the preferred model in CWHFTIO would 
require spherical shapes, which are very rare among their observed 
objects. In this section, we focus on tracing the cause of the dis¬ 
agreement between the results of the CWHFTIO analysis and the 
results of our own tests. 


3.1 A joint treatment of HI and molecular clouds? 

One likely culprit is the set of assumptions regarding the family 
of models CWHFTIO used to describe their data. The conclu¬ 
sion of CWHFTIO about the B — p relation is obtained through 
a careful Bayesian analysis. Despite the important conceptual 
and methodological advantages of Bayesian statistics in treat¬ 
ing diverse datasets, Bayesian analyses are necessarily model- 
dependent: they select the best parameter values in a specific 
family of models to treat the data at hand; but they do not con¬ 
vey information on how good a fit to the data the adopted family 
of models is as a whole. For example, a Bayesian analysis can 
always find a set of parameters of a normal distribution that best 
describe a specific dataset, even if the data are not distributed 


normally. However, that “best description” is actually a very poor 
representation of reality. 

A central assumption in the CWHFTIO family of models 
that could be affecting their findings on the value of the expo¬ 
nent n is that concerning low densities (HI data), namely, that 
the average magnetic field remains constant with density. This as¬ 
sumption is problematic for two reasons. First, observationally, it 
is not consistent with our empirical understanding of the behav¬ 
ior of magnetic fields with density for HI clouds. Several values 
of n have been proposed, but the value = 0 has not been one 
of them (e.g., see review by Vallee 1997). This has also been re¬ 
ported by Marchwinski et al. (2012), who explicitly contrasted 
their results to the model adopted by CWHFTIO at low densi¬ 
ties. Second, if one forces data on the B-p plane on a horizontal 
line at low densities and demands continuity of the B-p relation 
between low- and high-density datapoints, one necessarily sets a 
pivot point for the high-density part of the relation. If the value 
of the magnetic field at the “transition density” no (the density 
which separates the low-density and high-density branches of the 
model) is too low, the value of k, for the high-density part of the 
B-p relation will be forced to acquire greater values in order to 
accommodate the magnetic fields at the highest densities. 

We devise a “goodness of fit” test to quantify whether indeed 
the assumption that the average magnetic field in HI clouds is 
independent of density is consistent with the data, and whether 
this assumption can affect the global model fit. 

We start with the set of observed number densities for the 
objects used in the CWHFTIO analysis. For each number density 
rii, we randomly select a magnetic field according to the best 
generalized model by CWHFTIO. First we calculate the Bmax 
value appropriate for ni through. 


-Bmax ('?T') 


Bo, n < no 

So(^)“, n>no 


( 2 ) 


where Bo = 10 /rC, no = 300 cm“^, a = 0.65. We then select 
a total magnetic field Bi of the cloud from a uniform distribu¬ 
tion with boundaries between / • Bmax and Bmax, as prescribed 
by CWHFTIO. To account for various orientations we multiplied 
each Bi by cos(0) randomly drawn from a uniform distribution be¬ 
tween — I and 1, to obtain the line-of-sight magnetic field Blqs,!- 
We then “observed” this Blos,! W drawing a random value from 
a normal distribution with mean equal to Blos,! and standard de¬ 
viation equal to the actual observational uncertainty (Tb recorded 
by CWHFTIO for the cloud with density ni. 

In this way, we produce “mock” magnetic field observations 
at the same densities as the CWHFTIO dataset that are consis¬ 
tent with the CWHFTIO model for the magnetic field. We then 
compare the mock data with the actual magnetic field measure¬ 
ments through a Kolmogorov-Smirnov test. We perform the test 
for (i) the entire dataset, (ii) just HI observations, and (iii) just 
molecular observations. When we treat the entire dataset (HI and 
molecular observations), the K-S p-value for the a = 0.65 case 
was 19.7% - the model is acceptable. However, when we treat the 
HI and molecular data separately, the picture changes: for the 
molecular observations alone (where B is dependent on p accord¬ 
ing to the CWHFTIO model), the p-value for a = 0.65 is only 
5.2% - only marginally consistent with the data. In the case of 
HI data, the K-S test p-value is only 0.35%: B oc p*^ is not a good 
description of the low-density data. 

Figure [7]shows the cumulative distribution functions (CDFs) 
of the measured line-of-sight value of the magnetic field, for the 
different datasets, and demonstrates what happens when we com¬ 
bine HI and molecular data: both HI and molecular datasets are 
in poor agreement with the model, but in opposite directions: one 
produces too few low values of B (data CDF starts out below the 
model), and the other produces too many (data CDF starts out 
above the model). Adding the two datasets together moves the 
data CDF closer to the model. 
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Figure 7. Cumulative distributions of the actual observations (dashed blue lines) and mock-observed values drawn from the CWHTFIO 
“best-fit” model (green solid lines) of each of the following cases: entire dataset (left), just molecular observations (middle), and just HI 
observations (right). 


We have therefore shown quantitatively in this section that 
treating HI and molecular clouds jointly can affect the accept¬ 
ability of a value of k for the molecular branch of the model. For 
this reason, for the remainder of this work, we treat the molecular 
data by themselves. 


3.2 An uncertainty of only a factor of 2 in volume 
densities? 

An additional potentially problematic issue in the CWHFTIO 
analysis is the estimate of uncertainties in the density of each 
cloud/core. In CWHFTIO the issue is stated as follows: “Based 
on our experience in making such estimates, we choose to trust 
the inferred value within a factor of two, although the actual 
degree of uncertainty is not precisely known” (Cruether et al. 
2010). However, the uncertainty in the data is an important fac¬ 
tor in any statistical analysis. Underestimated uncertainties can 
make a dataset appear much more constraining than it really is. 
For this reason, we try to assess the actual degree of uncertainty 
for as many objects as possible by conducting a literature search 
of density estimates for each object and examining the spread 
between different estimates. 

In approximately half of the total number of molecular data- 
points given in Crutcher et al. (2010), volume densities are taken 
from surveys different from those from which the magnetic field 
data are taken. Even if the centre coordinates and angular resolu¬ 
tion between the Zeeman survey and the one probing other physi¬ 
cal parameters are the same, volume densities are usually derived 
from CO and CS measurements (see appendices in Crutcher 1999, 
and Falgarone et al. 2008). This introduces a bias since these 
molecules provide information for regions different from those 
where Zeeman emission occurs. Additional chemical assumptions, 
as for example the N(H 2 )/N(CO) ratio and the depletion of the 
species onto dust grains, can further complicate the picture (see 
Tassis et al. 2012; 2014). 

An example of the spread of volume-density estimates using 
different tracers is DR210H: Vallee & Fiege (2006) used ^^CO 


observations to determine the volume density of the core, which 
they reported to be n(H 2 ) = 3.4 x 10^ cm“^. Mangum et al 
(1991) used C^^O observations to derive n(H 2 ) > 6 x 10® cm“® 
for the same object. The difference is over an order of magnitude, 
much greater that the factor of 2 used by CWHFTIO. Further 
uncertainties are introduced through the morphology of the cores. 
In absence of a better practice, column densities are converted 
to volume densities by assuming spherical geometry. However, 
no real information exists for the line-of-sight size of the core. 
It is clear from the above that the uncertainties in the volume 
densities can be much greater than the factor of two adopted by 
CWHFTIO. 

Studying the cited literature, we find that volume-density 
values reported in CWHFTIO are not always consistent with 
the references to which they are attributed, presumably as a re¬ 
sult of updated data analysis (cores NGC2024, L1457S, L1457Sn, 
L1495(6), L1534, L1544). In any case, the changes in the adopted 
volume densities are indicative of the uncertainties introduced by 
various choices in the data analysis. 

Even if discrepancies greater than a factor of two are only 
present in a small number of cases, since the total number of 
points is small, the resulting effect on the statistical analysis of 
the dataset can be significant. We find the greatest variations 
for the more evolved cores; measurements of volume densities in 
other literature sources are generally greater than the ones re¬ 
ported in CWHFTIO. Table ^ summarizes different volume den¬ 
sity measurements for the fraction of CWHFTIO objects for which 
a literature search yielded additional density estimates. Figure 0 
demonstrates the effect of density uncertainties on the slope of the 
B - p relation in a simple way, focusing on datapoints with signif¬ 
icant (> 3cr) measurements of the magnetic-field strength rather 
than upper limits. Objects for which multiple measurements of 
the volume density are available in the literature are shown in 
color. If the lowest available measurements of the volume density 
are used, the value of (apparent slope derived through log-log 
regression on the B - p plane using detection-only data) is kw 
= 0.58 (solid line). These points are consistent with the volume- 
density values adopted by CWHFTIO. If the highest available 
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Figure 8. Two least-squares fits to detections only, when upper and lower limits of volume density measurements (curves labeled 
by /Cul and respectively) are adopted from the literature. Even when Zeeman measurements are accurate enough, volume density 
uncertainties can produce ambiguous results. Objects for which multiple measurements are available, and references, are given in Table 

m 


measurements of the volume density are used, then = 0.46 
(dashed line). 

We can also quantify the statistical effect of the updated 
uncertainties, by performing a Monte-Carlo simulation of the ap¬ 
parent slope hx derived through a naive statistical analysis of mock 
observations drawn from the CWHFTIO model with updated un¬ 
certainties in the volume density where these are available (and 
a factor of 2 when no additional information is available). Our 
analysis is similar to that of Crutcher (1999): a single power law 
is fitted through regression to the detections-only data (as de¬ 
fined in Crutcher 1999: measured line-of-sight \B\ in units of its 
observational uncertainty greater or equal than 2.5). We repeat 
the “experiment” 10^ times. The distribution of k values obtained 
in this manner is shown in Figure O A Gaussian fit to the dis¬ 
tribution shown yields a mean of 0.59 and a spread of 0.07: the 
detections-only slope is nonconstraining, and can be anywhere 
from 0.4 to ~ 0.8 within 2cr (assuming the CWHFTIO model 
is correct). 

We conclude that there are several cases for which the true 
uncertainty in the volume density is much greater than the factor 
of 2 adopted by CWHFTIO, and that the effect of changing the 
adopted volume-density measurements within the range of esti¬ 
mates found in the literature can have a very significant effect on 
the B-p relation of at least detections-only data. 

In order to quantify the effect of the large volume-density 
uncertainties in the treatment of the entire datasets, we re¬ 
peat the goodness-of-fit test of ^3.11 for molecular objects only, 
properly accounting for uncertainties in density: for objects for 
which multiple volume-density estimates exist in the literature, 
we choose a density uniformly distributed in the available range. 
For other objects, we assign a density within a factor of two of 
the CWHFTIO value. We calculate the cumulative distribution of 
simulated (mock-observed) magnetic-field values and we check its 
consistency with the cumulative distribution of actual magnetic- 
field strength observations through a K-S test. The hypothesis 
that the two distributions are the same is rejected at the 0.3% 
level. 



Figure 9. Distribution of apparent slope k, derived through log- 
log regression on the B-p plane of detection-only data for the 
optimal model of CWHFTIO, using updated volume-density un¬ 
certainties. 

The reason for this discrepancy can be seen in Figure HOI 
which overplots data and model predictions on the B-p plane. 
Using the same algorithm described above, we calculate, at each 
density bin, the median simulated \Bz\ (solid line) as well as the 
1^2:1 limits that contain Icr (68%, dashed lines) and 2cr (95%, 
dotted lines) of our simulated observations. On the same plot, we 
overplot with red dots the measurements that come from molec- 
































































10 Tritsis et al. 


Table 1. Density measurements of various objects in CWHFTIO. 
Values from: Crutcher et al. 1999 (1); Mookerjea et al 2012 (2); 
Troland & Crutcher 2008 (3); Hirota et al. 1997 (4); Motte et al. 
2007 (5); Girart et al 2013 (6); Sarma et al 2013 (7); Saito et al. 
2007 (8). 


Core 

CWHFTIO value 
(cm-3) 

Other Measurements 
(cm-®) 

Ref. 

NGC 2024 

1.3 X 10^ 

2 X 10^ 

(1) 

S88B 

1.3 X 10^ 

3 X 10^ 

(8) 



3.4 X 10^ 

(7) 

DR210HI 

3.4 X 10® 

2 X 10® 

(1) 



1.36 X 10® 

(2) 



2 X 10® 

(5) 



2 X lO'^ 

(6) 

DR2IOH2 

3.4 X 10® 

4 X 10® 

(1) 



1.36 X 10® 

(2) 



2 X 10® 

(5) 



2 X lO'^ 

(6) 

LI457S 

1.3 X 10-^ 

3 X 10® 

(3) 

LI457Sn 

1.7 X 10“^ 

4 X 10® 

(3) 

LI534 

5.4 X 10® 

I X 10^ 

(3) 

LI544 

3.2 X 10® 

1.3 X 10^ 

(3) 

L1595(6) 

4.1 X 10® 

8.2 X 10® 

(3) 

B1 

2 X 10^ 

3.2 X 10® 

(4) 


ular tracers in CWHFTIO. The cause of the strain registered by 
the K-S test between the CWHFTIO model and the data can be 
clearly identified in this plot. While reasonable fractions of points 
are within the expected Icr and 2(7 limits at various densities, they 
are not symmetric about the expected median: a larger fraction 
of points lies systematically below the median than above. 


3.3 A flat distribution of magnetic fleld strengths? 

A final potentially problematic assumption in CWHFTIO is that 
of the shape of the probability distribution of magnetic-field val¬ 
ues. CWHFTIO adopt a uniform distribution of H-values, be¬ 
tween some minimum and some maximum value. There are two 
issues arising from this choice. The first one is statistical. The uni¬ 
form distribution is a convenient tool and significantly simplifies 
statistical analyses when it can be used without loss of generality 
or in absence of further information. However, it is always the 
most restrictive option among all frequently used bi-parametric 
distributions of the same variance cr^. The reason is that, for a 
uniform distribution of finite width, the probability to have in¬ 
trinsic (i.e., free of observational error) values of the quantity of 



Figure 10. Median (solid line), Icr (dashed lines) and 2cr (dot¬ 
ted lines) limits of the simulated \Bz\ at each density according 
to the CWHTFIO model, overplotted with molecular data from 
CWHTFIO. 


interest at a distance greater than y/Sa from the mean is exactly 
zero. This means that this type of distribution cannot adequately 
describe peaked distributions with tails. If a collection of data- 
points that shows a preferred value (a peak) as well as outliers 
(tails) is forced onto a uniform distribution of variable width, the 
best fit that will be obtained by likelihood analysis is a uniform 
distribution that completely misses the peak and stretches out 
all the way to the farthest outlier. In terms of the B - p scal¬ 
ing problem: in order to accommodate, say, a single abnormally 
high value of the B-field at very high densities, the uniform B- 
distribution has to stretch its maximum to very high values, and 
this could again result in a steeper scaling than otherwise war¬ 
ranted by the data. This effect of the uniform-distribution choice 
is somewhat moderated by the treatment of observational uncer¬ 
tainties (making the probability of observation of an outlier finite 
due to observational error), but it is not clear that this is enough 
to eliminate the bias that could potentially be introduced. 

The second problem with the choice of a uniform distribution 
of H-values is conceptual. CWHFTIO assign a physical interpre¬ 
tation to the width of the best-fit uniform distribution they derive 
from their data: that, since the best-fit distribution extends uni¬ 
formly from almost zero to a maximum value, this implies that 
there is no preferred value of the magnetic field in objects of a 
specific density, and that this is additional evidence for the dy¬ 
namical insignificance of magnetic fields, compounded with their 
preferred 2/3 slope of the B - p scaling. However, if there were a 
preferred value of the magnetic field in the interstellar medium, 
but there were also outliers, the CWHFTIO uniform distribution 
would miss it and, consequently, would lead to a misleading phys¬ 
ical interpretation of the available data. 

The first question we ask in addressing these issues is 
whether a different family of distributions which is capable of de¬ 
scribing peaked distributions with tails may be an equally good or 
better fit to the magnetic-field data. Such a family of distributions 
is the lognormal. 


p{B) 


1 

-^— exp 


(InB - In Hm,a)^ 

2al 


(3) 
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which represents a Gaussian distribution of In B centered at 
lnBm,a and with spread cro- We will identify optimal parame¬ 
ters for a lognormal (instead of a uniform) distribution of in¬ 
trinsic magnetic-field strengths, which is consistent with the 
observed \Bz\ values in a narrow density range (from 2 x 10^ to 
4 X 10^ cm“^). We have chosen that particular density range be¬ 
cause it is well populated (27 observed objects) so that it provides 
enough statistics to give a good sense of the underlying distribu¬ 
tion, and narrow enough so that we can ignore any evolution with 
density between points. 

In order to determine parameters of this distribution, we 
scan the two-dimensional parameter space of (Sm,a, <^ 0 ), and for 
each pair we simulate a distribution of observed | Bz | with random 
observation directions (uniform cos 0) and including Gaussian ob¬ 
servational errors identical to the ones quoted in GWHTFIO. The 
resulting distribution is then compared to the distribution of ob¬ 
served \Bz\. Since the purpose of this work is not to repeat the 
sophisticated statistical analysis of GWHTFIO but simply to as¬ 
sess the effect of relaxing the potentially problematic assump¬ 
tions, we have not formally obtained a best-fit distribution; in¬ 
stead, we selected the set of parameters that minimizes the K-S 
statistic (yielding, in our case, a K-S p-value of 80% of the two 
datasets to be drawn from the same distribution). This means 
that a proper fit (for example, through a maximum-likelihood 
analysis) may yield a slightly different set of optimal parame¬ 
ters, which, if anything, will be an even better fit to the data. 
The optimal parameters we have identified for the lognormal in 
this narrow density range are \n(Bm,a/= 2.57 and ao = 0.3, 
respectively. We assign this distribution to a number density of 
ria = 2.61 X 10^cm“^, which is the average of the GWHTFIO 
quoted densities for objects in the density range we have consid¬ 
ered. In the analysis that follows we assume that, for all densities 
above 300 cm“^, the distribution p{B) remains lognormal with 
the same ao at all densities, while Bm scales as 

-Hm('^) = Hm,a('^a) ( ^ • (4) 

\naj 

Figure ITT] shows how the observed \Bz\ values in the density 
range we have considered compare to those expected from the 
optimal lognormal B distribution and the GWHTFIO uniform 
B distribution (assuming in both cases uniform viewing angles 
cos 6 and observational uncertainties as quoted in GWHTFIO). It 
is clear that the lognormal is a better match for the qualitative 
behavior of the observed data. 

Figure fl^ shows the corresponding probability density func¬ 
tions for the intrinsic magnetic-field strength B (instead of the 
observed line-of-sight strength of the magnetic field, \Bz\) at the 
number density ria- The plot demonstrates how the GWHFTIO 
conclusion, that there appears to be no preferred value of the 
magnetic field at a specific density in the interstellar medium, 
because a wide uniform distribution is preferred over a much nar¬ 
rower one, is a direct result of the constraining nature of a uni¬ 
form distribution: the optimal lognormal, which we have shown 
is a better qualitative description of the observed data, does show 
a significant peak, corresponding to a preferred value at a given 
density; however, significant tails exist. In contrast, the optimal 
uniform distribution misses this peak, and extends to high and 
low values in order to accommodate the tail values, giving the 
false impression that the spread in intrinsic B needed to explain 
the observed Bz is much greater than it is in reality. 


4 RECONCILING THE B-p RELATION 
WITH CORE SHAPES. 

Having identified three potentially problematic assumptions in 
the GWHFTIO analysis, we now re-evaluate the information their 
{Bz, n) datapoints convey regarding the slope of the B-p relation 



|BJ (nG) 

Figure 11. Distribution of line-of-sight magnetic field values Bz 
in the density range 2 — 4 x I0^cm“^. Blue histogram (solid 
line): all observed objects (molecular data only). Green histogram 
(dashed line): optimal lognormal distribution. Red histogram 
(dashed-dotted line): GWHTFIO model (uniform distribution). 


and we reconcile the latter with the lack of evidence for a prefer¬ 
ence for spherical geometry in the GWHFTIO cores. We do so by: 
(a) treating molecular (high-density) data on their own, since this 
is where the debated scaling (B oc vs B oc p^/^) arises; (b) 
using updated uncertainties in the volume density where these 
are available; and (c) using a lognormal distribution of intrinsic 
magnetic-field strengths. 

We repeat our goodness-of-fit tests of ^34] using the optimal 
lognormal distribution described in 3331 with its mean scaled 
with n with a slope a equal to either 0.65 (the GWHFTIO pre¬ 
ferred value) or 0.5 (the historically preferred value due to the¬ 
oretical expectations from magnetically-controlled gravitational 
contraction and due to the empirical results from fitting a power 
law to detection-only data). Figures [13] and im show the \Bz\ - 
n plane with lines corresponding to the median (solid), Icr (long 
dashes), and 2a (short dashes) expected limits of mock obser¬ 
vations drawn from the lognormal model, with slope a equal to 
0.65 and 0.5, respectively. A K-S test between the observed and 
simulated distributions of \Bz\ values in the entire density range 
of molecular datapoints, as in 3331 returns a p-value of 7.5% for 
observed and simulated data to be drawn from the same distri¬ 
bution, for a scaling slope of 0.65; it returns a p-value of 15.2% 
for a scaling slope of 0.5. 

The agreement between data and model has improved sig¬ 
nificantly even for a scaling slope of 0.65 (from 0.3% to 7.5%) 
by switching from a uniform distribution of B-values to a lognor¬ 
mal. This is achieved without changing the number of modeling 
parameters; both uniform and lognormal are bi-parametric fam¬ 
ilies of distributions. However, a slope of o; = 0.5 is preferred 
over a = 0.65, although the two cannot be distinguished at a 
statistically significant level. 
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Figure 12. Distribution of intrinsic magnetic-field strength B at 
n = 2.61 X 10^cm“^. Green (dashed) line: optimal lognormal; 
Blue (solid) line: CWHTFIO model (uniform distribution). 
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Figure 14. Median (solid line), Icr (long-dashed lines) and 2a 
(short-dashed lines) limits of the simulated \Bz\ at each density 
obtained from the optimal lognormal p{B) and a = 0.5, overplot¬ 
ted with molecular data from CWHTFIO. 



Figure 13. Median (solid line), Icr (long-dashed lines) and 2cr 
(short-dashed lines) limits of the simulated \Bz\ at each density 
obtained from the optimal lognormal p(B) and a = 0.65, over¬ 
plotted with molecular data from CWHTFIO. 


5 SUMMARY 

In this paper, we discussed ways to assess the relation between 
the magnetic-field strength and the gas density in the interstellar 
medium. 

We reviewed the connection between the exponent k in the 
relation B oc p'^ and the geometry of a cloud. We showed that 
specific combinations of cloud geometry and magnetic-field ori¬ 


entation result in different B - p relations, with B oc being 
unique to the spherical-cloud (or core) geometry. 

In light of this result, we sought to verify whether the claim 
of H oc (at high densities) through the statistical analysis 

of a large sample of density and magnetic-field measurements by 
CWHFTIO is consistent with the geometry of objects in their 
sample. To this end, we have studied emission maps in all 27 ob¬ 
jects in the CWHFTIO sample for which such data were available. 
We only found aspect ratios consistent with a spherical geometry 
in 4 of them; the distribution of aspect ratios does not show any 
preference toward unity, which would be the signature of spherical 
shapes. We thus concluded that there is no evidence of preferen¬ 
tially spherical objects in the CWHFTIO, a result inconsistent 
with a B oc p^/^ relation. 

We then investigated the possibility that this disagreement 
could be caused by simplifying assumptions in the CWHFTIO 
statistical analysis. 

We first tested the effect of using a joint model for low- 
and high-density data and forcing the low-density data to a con¬ 
stant magnetic-field strength, independent of density. We used 
Monte-Carlo simulations to produce mock observations from the 
best CWHFTIO model with identical observational uncertainties 
as the ones quoted in CWHFTIO, and we used a Kolmogorov- 
Sminrov test to determine whether the mock observations were 
consistent with the data. We found that: 

• When treating the two branches (HI and molecular) of the 
CWHFTIO model separately, the low-density B oc p° branch is 
inconsistent with the data (p-value = 0.35%), while the high- 
density branch is marginally consistent (p-value = 5.2%). 

• These two branches deviate from the data in opposite direc¬ 
tions, and, as a result, when we treat the combined dataset, we 
find artificially improved consistency (p-value = 19.7%). 

We therefore conclude that the finding of CWHFTIO, that 
a model with B oc p° at low densities and B oc p^/^ at high 
densities is the preferred description of observations of densities 
and magnetic fields, is an artifact of their combined treatment of 
low- and high-density datapoints. 
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We also checked whether the uncertainty in volume densities 
of a factor of 2 adopted by CWHFTIO is a good estimate, by trac¬ 
ing literature sources and comparing different density estimates 
for the same objects when these were available. We found that for 
several objects the actual uncertainties, as reflected in the spread 
of estimates in the literature, are much greater. When repeating 
our consistency analysis between the high-density branch of the 
CWHFTIO model and the data using updated volume-density 
uncertainties, the agreement between model and data worsens, 
with the p-value dropping to 0.3%. The reason for this is that the 
additional density values available in the literature (especially at 
the highest-density objects) tend to be greater, instead of being 
symmetrically distributed about the value adopted in CWHFTIO: 
more recent estimates have generally produced upward correc¬ 
tions in volume densities. 

We investigated whether a lognormal distribution for p{B) 
would yield better agreement with the data than the CWHFTIO 
uniform distribution. We found that it does. 

Finally, we relaxed the three problematic assumptions, by (a) 
treating molecular observations on their own, (b) using updated 
volume-density uncertainties, and (c) using the optimal longormal 
to model the distribution of magnetic-field strengths. We repeated 
our goodness-of-flt K-S tests for the lognormal models, scaled with 
density with slopes of 0.65 (the optimal CWHFTIO slope) and 0.5. 
We found that the K-S test accepts both models. However there 
is a preference for B oc (by a factor of 2 in the p-value of the 
K-S test), which is also preferred by the independent analysis of 
cloud shapes. This result is in agreement with predictions of the 
ambipolar-diffusion theory of star formation. 

ACKNOWLEDGEMENTS 

KT acknowledges support by FP7 through Marie Curie Career 
Integration Grant PCIG- GA-2011-293531 SFOnset. AT, GVP, 
and KT acknowledge partial support from the EU FP7 Grant 
PIRSES-GA-2012-31578 EuroCal. GVP, KT, and TM acknowl¬ 
edge support by the RoboPol project, which is implemented under 
the ARISTEIA Action of the OPERATIONAL PROGRAMME 
EDUCATION AND LIFELONG LEARNING and is co-funded 
by the European Social Fund (ESF) and Greek National Re¬ 
sources. TM is grateful for the hospitality of the Department of 
Physics of the University of Crete during part of this work and 
writing of this paper. AT, GP, and KT thank the Astronomy De¬ 
partment of Caltech, for its hospitality during various stages of 
this work. We have made ample use of the public Herschel Science 
Archive data products. We thank P. Goldsmith for providing the 
CO data of the Taurus Molecular Cloud and for comments on the 
manuscript. We are grateful to T. Hartquist for careful reading 
of the paper and insightful comments. 

REFERENCES 

Basu, S., Ciolek, G. E., 2004, ApJ, 607, 39 

Chapman, N. L., Davidson, J. A., Goldsmith, P. F., et ah, 2013, 
ApJ, 770, 151 

Ciolek, G.E., Basu, S., 2006, ApJ, 652, 442 
Crutcher, R. M., 1999, ApJ, 520, 706 

Crutcher, R. M., Wandelt, B., Heiles, C., Falgarone, E., Troland, 
T. H., 2010, ApJ, 725, 466 

Falgarone E., Troland T. H., Crutcher R. M., Paubert G., 2008, 
A&A, 487, 247 

Fiedler, R. A., Mouschovias, T. Ch., 1992, ApJ, 391, 199 
Fiedler, R. A., Mouschovias, T. Ch., 1993, ApJ, 415, 680 
Gammie, C. F., Lin, Y. T., Stone, J. M., Ostriker, E. V., 2003, 
ApJ, 592, 203 

Girart, J. M., Rao, R., Marrone, D. P., 2006, Science, 313, 812 
Girart, J. M., Frau, P., Zhang, Q., Koch, P. M., Qiu, K., Tang, 
Y. W., Lai, S. R, Ho, P. T. R, 2013, ApJ, 772, 69 


Goldsmith P. F., Heyer M., Narayanan J., Snell R., Li D., Brunt 

C. , 2008, ApJ, 680, 428 

Goodwin, S. P., Ward-Thompson, D., Whitworth, A. P., 2002, 
MNRAS, 330, 769 

Hirota, T., Yamamoto, S., Mikami, H., Ohishi, M., 1997, ApJ, 
503, 717 

Jones, C. E., Basu, S., 2002, ApJ, 569, 280 

Jones, C. E., Basu, S., Dubinski, J., 2001, ApJ, 551, 387 

Kudoh, T., Basu, S., 2011, ApJ, 728, 123 

Li, H.-b., Yuen, K.H., Otto, F., Leung, P.K., Sridharan, T.K., 
Zhang, Q., Liu, H., Tang, Y.-W., Qiu, K., 2015, DOI:10.1038 
Naturel4291 

Li, P. S., Norman, M. L., Mac Low, M. M., Heitsch, F., 2004, 

ApJ, 605, 800 

Mac Low, M. M., Klessen, R., S., 2004, RvMP, 76, 125 
Mangum, J. G., Wootten, A., Mundy, L. G., 1991, ApJ, 378, 576 
Marchwinski, R. C., Pavel, M. D., Clemens, D. P., 2012, ApJ, 
755, 130 

Mestel, L., 1965, QJRAS, 6, 265 
Mestel, L., Spitzer, L. Jr., 1956, MNRAS, 116, 503 
Mookerjea, B., Hassel, G. E., Gerin, M., Giesen, T., Stutzki, J., 
Herbst, E., Black, J. H., Goldsmith, P. F., Menten, K. M., 
Krelowski, J., De Luca, M., Csengeri, T., Joblin, C., Kazmier- 
czak, M., Schmidt, M., Goicoechea, J. R., Cernicharo, J., 
2012, A&A, 546, 75 

Motte, F., Bontemps, S., Schilke, P., Schneider, N., Menten, K. 

M., Broguiere, D., 2007, A&A, 476, 1243 
Mouschovias, T. Ch., 1976a, ApJ, 206, 753 
Mouschovias, T. Ch., 1976b, ApJ, 207, 141 
Mouschovias, T. Ch., 1979, ApJ, 228, 475 

Mouschovias, T. Ch., 1987a, in Physical Processes in Interstellar 
Clouds, ed. G. Morflll & M. Scholer (Dordrecht: Reidel), 453 
Mouschovias, T. Ch., 1991a, ApJ, 373, 169 

Mouschovias, T. Ch., 1991b, in The Physics of Star Formation 
and Early Stellar Evolution, ed. C. J. Lada & N. D. Kylafls 
(Dordrecht: Kluwer), 61 

Mouschovias, T. Ch., 1991c, in The Physics of Star Formation 
and Early Stellar Evolution, ed. C. J. Lada & N. D. Kylafls 
(Dordrecht: Kluwer), 449 

Mouschovias, T. Ch., Ciolek, G., E., Morton, S., A., 2011, MN¬ 
RAS, 415, 1751 

Mouschovias, T. Ch., Morton, S., A., 1991, ApJ, 371, 296 
Mouschovias, T. Ch., Morton, S., A., 1992a, ApJ, 390, 144 
Mouschovias, T. Ch., Morton, S., A., 1992b, ApJ, 390, 166 
Mouschovias, T. Ch., Spitzer, L., Jr., 1976, ApJ, 210, 326 
Myers, P., C., Fuller, G., A., Goodman, A., A., Benson, P., J., 
1991, ApJ, 376, 561 
Ostriker, J., 1964, ApJ, 140, 1056 

Pilbratt, G. L., Riedinger, J. R., Passvogel, T., Crone, G., Doyle, 

D. , Gageur, U., Heras, A. M., Jewell, C., Metcalfe, L., Ott, S., 
Schmidt, M., 2010, A&A, 518, LI 

Ryden, B. S., 1996, ApJ, 471, 822 

Saito, H., Saito, M., Sunada, K., Yonekura, Y., 2007, ApJ, 659, 
459 

Sarma, A. P., Brogan, C. L., Bourke, T. L., Eftimova, M., Troland, 
T. H., 2013, ApJ, 767, 24 
Tassis, K., 2007, MNRAS, 379, L50 
Tassis, K., Mouschovias, T. Ch., 2007, ApJ, 660, 402 
Tassis K., Dowell C.D., Hildebrand R.H., Kirby L., Vaillancourt 
J.E., 2009, MNRAS, 399, 1681 

Tassis, K., Willacy, K., Yorke, H. W., Turner, N. J., 2012, ApJ, 
754, 6 

Tassis, K., Willacy, K., Yorke, H. W., Turner, N. J., 2014, MN¬ 
RAS, 445, L56 

Troland T. H., Crutcher R. M., 2008, ApJ, 680, 457 
Vallee, J. P., 1997, Fund. Cosm. Phys., 19, 1 
Vallee, J. R, Fiege, J. D., 2006, ApJ, 636, 332 
Verschuur, G. L., 1969, Nature, 223, 140 




